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ABSTRACT 

We present radiative transfer models of the circumstellar environment of classical T Tauri 
stars, concentrating on the formation of the Ha emission. The wide variety of line profiles 
seen in observations are indicative of both inflow and outflow, and we therefore employ a 
circumstellar structure that includes both magnetospheric accretion and a disc wind. We per- 
form systematic investigations of the model parameters for the wind and the magnetosphere 
to search for possible geometrical and physical conditions which lead to the types of profiles 
seen in observations. We find that the hybrid models can reproduce the wide range of profile 
types seen in observations, and that the most common profile types observed occupy a large 
volume of parameter space. Conversely, the most infrequently observed profile morphologies 
require a very specific set of models parameters. We find our model profiles are consistent with 
the canonical value of the mass-loss rate to mass-accretion rate ratio (fi = 0.1) found in earlier 
magneto-hydrodynamic calculations and observations, but the models with 0.05 < \i < 0.2 
are still in accord with observed Ha profiles. We investigate the wind contribution to the line 
profile as a function of model parameters, and examine the reliability of Ha as a mass ac- 
cretion diagnostic. Finally, we examine the Ha spectroscopic classification used by Reipurth 
et. al, and discuss the basic physical conditions that are required to reproduce the profiles in 
each classified type. 

Key words: stars: formation - circumstellar matter - radiative transfer - stars: pre-main- 
sequence 



1 INTRODUCTION 

T Tauri stars (TTS) a re young (> 3 x 10 6 yrs, 
Appenzelle r & Mundtl ll98St) low-mass objects, and are the 
progenitors of solar-type stars. Classical T Tauri stars (CTTS) 
exhibit strong Ha emission, and typically have spectral types of 
F-K. Some of the most active CTTS show emission in higher 
Balmer lines and metal lines (e.g., Ca II H and K). They also exhibit 
excess continuum flux in the ultraviolet (UV) and infrared (IR). 
Their spectral energy distribution and polarisation data suggest the 
presence of circumstellar discs, which play an impo rtant role in 
regulating dynamics of gas flows around CTTS (e.g. ICamenzind 

Many observational studies (e.g., Herbig 1962; Edwards et al. 
[ 1994 iKenvon et all 11994 iReipurth. Pedrosa & Laed Il996t 
Ulencar & Basril l200d) of CTTS line profiles have revealed 
evidence for both outward wind flows and inward accretion 
flows, characterised by the blue-shifted absorption features 
in Ha profiles and redshifted inverse P Cygni (IPC) pro- 
files respectively. Typical mass-lo ss rates of CTTS are about 
1CT 9 to lO- 7 M yr- 1 (e.g., iKuhil 1 1964 lEdwards et all 
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Il987t lHartigan. Edwards & Ghandourl Il995l) . and the mass- 
accretion rates are al so ab ou t 10~ 9 to 10~ 7 M l ^yr~ 1 (e.g. , 
i Kenvon & Hartmannl Il987l : iBertout. Basri & Bouvierl Il988t 
Gullbring et al. 1998) . Recent Ha spectro-astr ometric observa- 
tions bvlTakami. Bail ey & Chrvsostomod ^2003) provide indirect 
evidence for the presence of bipolar and monopolar outflows down 
to ~ 1 au scale (e.g. CS Cha and RU Lup). Similarly, ESO VLT 
observations using high-resolution (R = 50 000) two-dimensional 
spe ctra of edge-on C T TS (H H 30, HK Tau B, and HV Tau C) 
by I Appenzeller et alJ l2005t) show extended Ha emission in 
the direction perpendicular to the obscuring circumstellar disc, 
suggesting the presence of the bipolar outfl ows. On an e ven larger 
scale, HST observations of HH30 I Bu rrows et alJll99^) trace the 
jet to within < 30 au of the star. The jet has a cone shape with 
an opening angl e of 3° betwee n 70 and 700 au (Kon igl & Pudritzj 
2000). lAlencar & Basril ^2000) found about 80 per cent of their 
sample (30 CTTS) show blue-shifted absorption components in at 
least one of the Balmer lines and Na D. 

In the currently favoured model of accretion in CTTS, the ac- 
cretion disc is disrupted by the magnetosphere, which channels the 
gas from the disc onto the stellar surface ( e.g., Uchida & Shibata 
|l985 tlK6nigMl99ltlCollier Cameron & Campbell 1993tlShu et all 
1994). This picture is supported by recent measurements of strong 
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(~ 10 3 G) magnetic fiel ds in CTTS (e.g.. ljohns-Krull et alJll999t 
ISvmington et alj[2005bl) and by radiative transfer models which 
reprodu ce the gross ch aracteristi cs of observe d profiles for some 
CTTS {Muzerolle. Calvet & Hartmanr] [200 ll) . In particular the 
magnetospheric accretion (MA) model explains blue-ward asym- 
metric emission line profiles as resulting from the partial occulta- 
tion of the flow by the stellar photosphere and by the inner part of 
the accretion disc, while inverse P Cygni profiles in the MA model 
result from inflowing material at near free-fall velocities seen pro- 
jected against hotspots on the stellar surface. 

Despite these successes, the overwhelming observational evi- 
dence for outflows in CTTS suggests that the MA model is only one 
component of a complex circumstellar environment. Clearly, one 
must include the contribution of any wind/jet flow if one wishes 
to both accurately predict the mass-accretion rate and also de- 
termine the mass-loss rate of a CTTS via emission profile mod- 
elling. The first attempt in this direction was made bv lAlencar et alJ 
1 2005) who demonstrated that the observed Ha, H/3 and Na D lines 
of RW Aur are better reproduced by the radiative transfer model 
which included a collimated disc-wind arising from near the inner 
edge of the accretion disc. 

The main aim of this paper is to find a simple kinematic model 
which can reproduce the wide variety of the observed profiles, and 
to perform empirical studies of line formation in an attempt to place 
morphological classification schemes on a firmer physical footing. 
We will also discuss whether our model is consistent with some 
predictions made by recent (magneto-hydrodynamics) MHD stud- 
ies i.e. fi = Mwind/Macc « 0.1 (e.g. lKonigl & PudritzfeOOOh . 

In Section|2| the model assumptions and the basic model con- 
figurations are presented. We discuss the radiative transfer model 
used to compute the profiles in Section[3] and the results of model 
calculations are given in Section |4| We discuss our results in the 
context of Reipurth's classification scheme in Section [5] and our 
summary and conclusions are presented in Section|6| 



2 MODEL CONFIGURATION 

In order to understand how the different parts of the CTTS circum- 
stellar environment contribute to the formation of Ha, the model 
space is divided into four different regions: (1) a central continuum 
source, (2) the magnetospheric accretion flow, (3) the wind out- 
flow, and (4) the accretion disc. Fig.Qdepicts the relative location 
of the regions in the model space. The density is assumed to be ro- 
tationally symmetric around the z-axis. The innermost radius of the 
magnetosphere at the equatorial plane coincides with the the inner 
radius of the accretion disc. From the innermost part of the accre- 
tion disc, the gas falls freely, moving along the magnetic field onto 
the surface of the star. In the following subsections we describe the 
details of the model components. 



2.1 The continuum source 

We adopt stellar parameters of a typical classical T Tauri star for the 
central continuum source, i.e. radius (R,), mass (M*), and effec- 
tive temperature of the photosphere (T p h) are 2 R^, . 5 M^ , and 
4000 K respectively. The model atmosphere of Kurucz 1 1979) with 
r p h = 4000 K and log(?„ = 3.5 (cgs) defines the photospheric 
contribution to the continuum flux. The parameters are summarised 
in TableQ 

An additional continuum source is considered for models 



which include magnetospheric accretion: as the infalling gas ap- 
proaches the stellar surface, it decelerates in a strong shock, and is 
heated to ~ 10 6 K. The X-ray radiation produced in the shock will 
be absorbed by the ga s locally, and re-emitted a s optical and UV 
light iCalve t & Gullb rin3ll99a iGullbring etai]l200d) . This will 
create hot rings where the magnetic field intersects with the sur- 
face. We assume that the free-falling kinetic energy is thermalized 
in the radiating layer, and is re-emitted as blackbody radiation with 
a single temperature. With the parameters of the magnetosphere 
and the star given above (Table0, about 8 per cent of the surface 
is covered by the hot rings. In our model, the fraction of hot-ring 
coverage area is determine d only by the ge ometry of the magneto- 
sphere (Section ^, 2> , Gullbring et al . (2000) estimated the fractions 
of the hot-ring coverage area for three intermediate CTTS by using 
the observed UV and optical spectra. They found the range of frac- 
tion to be 0.3-3.0 per cent. This is smaller than the value adopted 
here; however, even if we use 1 per cent hot-ring area coverage 
model (by changing the width of the magnetosphere), we find the 
line sh apes of Ha do not change significantly (see also Fig. 5 of 
iMuzerolle et all200lh . 

If the mass-accretion rate is 10~ 7 Mq yr _1 , the ratio of the 
accretion luminosity to the photospheric luminosity is about 0.5, 
and the corresponding temperature of the hot rings is about 6400 K. 
The luminosity ratio and the temperature depend on an adopted 
mass-accretion rate. The continuum emission from the hot rings is 
taken into account when computing the line profiles. 

2.2 The magnetosphere 

We ad opt the MA flow model of [Hartmann^few ett & Calvel 
Jl994l). as done by IMuzerolle et alJ 1200 ll) and by 
Svmington, Harries & Kurosawa (2005a), in which the gas 
accretion on to the stellar surface from the innermost part of the 
accretion disc occurs through a dipolar stellar magnetic field. 
The magnetic field is assumed to be so strong that the gas flow 
does not affect the underlying magnetic field itself. As shown in 
Fig.Q(region 2), the innermost radius (i? m of the magnetosphere 
at the equatorial plane (z = 0) is assigned to be the same as 
the inner radius (Rdi) of the accretion disc where the flow is 
truncated. In our models, i? m i and the outer radius (R mo ) of the 
magnetosphere (at the equatorial plane) are set to be 2.2 R, and 
3.0 R, respectively. The geometry of the magnetic field/stream 
lines is fixed for all calculations. We note that this magnetospheric 
geome try is identical to the "small/wide" model of Muzerolle et al. 
I2001h . 

The magnetic field and the gas stream lines are assumed to 
have the following simple form 

r = _R m sin 2 8 (1) 

(see iGhosh. Pethick & Lamblll977h where r, and 8 are the radial 
and the polar component of a field point position vector in the ac- 
cretion stream. R m is the radial distance to the field line at the equa- 
torial plane (8 = 7r/2), and its value is restricted between f? m i and 
Rmo (Fig.0. Using the field geometry above and conservation of 
energy, the velocity and t he density of the accret ing gas along the 
stream line are found as in lHartmann et aljll994l) . 

The te mperature structure of the magnetosphere used by 
iHartmann et alJ 1 19941) is adopted here. They computed the tem- 
perature, assuming a volumetric heating rate which is proportional 
to r~ 3 , by solving the ene rgy balance of the r adiative cooling 
rate oflHartmann et alJ ll982h and the heating rate l Hartman n et alJ 
ll994t) . IMartinl ll996l) presented a self-consistent determination 
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Figure 1. Basic model configuration of the disc-wind-magnetosphere hybrid model. The system consist of four components: (1) the continuum source located 
at the origin of the cartesian coordinates (x, y, z) - the j/-axis is into the paper, (2) the magnetospheric accretion flow, (3) the disc wind, and (4) the accretion 
disc. The disc wind originates from the disc surface between Wi = R w i and Wi = ii wo where Wi is the distance from the z axis on the equatorial plane. The 
wind source points (S), from which the stream lines diverge, are placed at distance d above and below the star. The degree of wind collimation is controlled 
by changing the value of d. 



Table 1. Summar y of the reference class ical T Tauri star model parameters. Note that a typical mass-accretion rate of CTTS is ~ 3 X 10 Mq yr (e.g. 
IGullbring et al]|l998tlCa?vet et alj|2004 ; however, a slightly higher mass-accretion rate is adopted here to demonstrate the effect of line-broadening more 
clearly (c.f. Section l4**TV 



Parameters R, M» T ph i? m i i? mo M a cc Af wind R di R do 
[Rp] [Mq] [K] [Bj [R*] [M yr-!] [M 0yr - 1 ] [R,] [au] 

2.0 0.5 4000 2.2 3.0 10" 7 10" 8 2.2 100 



of the thermal structure of the inflowing gas along the dipole 
magnetic field (equation by solving the heat equation cou- 
pled to the rate equations for hydrogen. He found that main heat 
source is adiabatic compression due to the converging nature of 
the flow, and the major contributors to the cooling process are 
bremsstrahlung radiation and line emission from Ca II and Mg II 
tons. iMuzerolle. Calvet. & Hartmannl 1 1998) found that the line 
profile model s computed according to the temperature structure of 
Martin Il99r3l do not agree with observations, unlike profiles based 
on the (less self-consistent) Hartmann et al. 1 1994) temperature dis- 
tribution. It is clear that the temperature structure of the magneto- 
sphere is still a large source of uncertainty in the acc retion model; 
here we adopted the model of H artmann et aljfT994l) . 



2.3 The disc wind 

The magneto-centrifugal wind paradigm, first proposed by 
Blandf ord & Pavnd (l 982 ) , has been often used to model the large- 
scale wind structure of T Tauri stars, or t h e observed optic al jets 
(e.g. HH 30 jet by iBurrows et ail Il996t |Rav et alj Il99fj|) . The 
launching of the wind from a Keplerian disc is typically done 
by treating the equatorial plane of the disc as a mass-inj ecting 
boundary condition (e.g., |Shu_etalJLl994t Q jsjyjig ova et a l. 1995; 
lOuved & Pudritzll 19971 iKrasnopolskv et alj|2003l) . Depending on 
the location of the open magnetic fields anchored to the disc, two 
different types of winds are produced. If the field is constrained 
to be near the c o-rotation radius of the stellar magnetosphere, an 
"X-wind' ' IShu et alll99"4T) is produced. If the open field lines are 
located in a wider area o f the disc, a "disc-win d" similar to t hat of 
iKonigl & PudritJ fcOOd) is produced (see also lKrasnopolskv et alJ 
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2003). Rec ent reviews on th e jet/wind-disc connecti on can be found 
in lKonigl & Pudridi20obT) and lPudritz & Baneriedl2005l) . Clearly 
there are several alternative outflow scenarios. Here we adopt a sim- 
ple kinematical model, based on the disc-wind paradigm, which 
broadly represents the results of MHD simulations, 

iKnigge. Woods & Drew! 1 19951) introduced the "split- 
monopole" kinematic disc-wind model in their studies of the UV 
resonance lines formed in the winds of cataclysmic variable stars. 
Their formalism provides a simple parameterisation of a disc wind 
that has similar properties to those found by MHD modelling. 
In this model, the outflow arises from the surface of the rotating 
accretion disc, and has a biconical geometry. The specific angular 
momentum is assumed to be conserved along a stream line, and 
the poloidal velocity component is assumed to be simply a radial 
from "sources" vertically displaced from the central star. Here we 
briefl y describe the disc-wind model, an d read ers are referred to 
IKnigge et all < 19951) and lLong & Kniggd <2002h for details. 

The four basic parameters of the model are: (1) the mass-loss 
rate, (2) the degree of the wind collimation, (3) the velocity gradi- 
ent, and (4) the wind temperature. The basic configuration of the 
disc-wind model is shown in Fig.Q The disc wind originates from 
the disc surface, but the "source" point (5*), from which the stream 
lines diverge, are placed at distance d above and below the centre 
of the star. The angle of the mass-loss launching from the disc is 
controlled by changing the value of d. The mass-loss launching oc- 
curs between R w i and R wo where the former is set to be equal to 
the outer r adius of the magnetosphere (i? mo ) and the latter is set to 
1 au, as in lKrasnopolskv et alJl2003l) . 

The local mass-loss rate per unit area (m) is assumed to be 
proportional to the mid-plane temperature of the disc, and is a func- 
tion of the cylindrical radius w = (x 2 



2\l/2 

y ) ,i-e. 



rh(w) oc [T(w)] 4 



(2) 



The mid-plane temperature of the disc is assumed to be expressed 
as a power-law in w; thus, T oc w q . Using this in the relation above, 
one finds 



m (w) oc w p 



(3) 



where p — 4a x ?. The index of the mid-plane temperature 
power law is adopted from the dust radiative transfer model of 
Whit ney" et alJ 120031) who found the innermost part of the accre- 
tion disc has q = —1.15. I n order to be consis t ent w ith the col- 
limated disc-wind model of Krasnooolskv et al. (2003) who used 
p — —7/2, the value of a is set to 0.76. The constant of propor- 
tionality in equation|3|is found by integrating rh from R w i to R wo , 
and normalising the value to the total mass-loss rate M W md- 

The azimuthal/rotational component of the wind velocity 
i>0 (w, z) is computed from the Keplerian rotational velocity at the 
emerging point of the stream line i.e. v<j, (wt, 0) = (GM, /wi) 1 ' 2 
where Wi is the distance from the rotational axis (z) to the emerging 
point on the disc, and by assuming the conservation of the specific 
angular momentum along a stream line: 



v<t,(w, z) =V<t,(Wi,0) (— ) 
V w I 



(4) 



Bas ed on the canonical velo city law of hot stellar winds 
(c.f. lCastor. Abbott & Kfeirll975l) . the poloidal component of the 
wind velocity (v p ) parameterised as: 



(5) 



V p (Wi,l) = C s (Wi) + [fv csc — C s (Wi)] ( 1 — — ^ j 

where c s , /, and I are the sound speed at the wind launching point 



on the disc, the constant scale factor of the asymptotic terminal 
velocity to the local escape velocity (from the wind emerging point 
on the disc), and the distance from the disc surface along stream 
lines respectively. R 3 i s the wind scale length , and its value is set 
to 10 i?wi by following lLong & Kniggd 120021) . 

Assuming mass-flux conservation and using the velocity field 
defined above, the disc wind density as a function of Wi and I can 
be written as 



p(wi,l) 



d 



(6) 



v p (wi,l) |cos 5| \D(iVi,l) cos 5 

where D and S are the distance from the source point (S) to a 
point along the stream line and the angle between the stream line 
and the disc normal, respectively. Fig.|2|shows the density and the 
velocity components along the mid stream line, passing through 
Wi — (R w i + i?wo) /2, on the disc plane (z = 0) for different 
values of the wind acceleration parameter f3. 

2.4 The accretion disc 

We adopt a simple analyt ical accretion disc model, the q-disc 
'stand ard model' IShakura & Sunvae^ll973l : lFrank. King & Raind 
2002 ) with the inner radius fixed at the inner radius of the magneto- 
sphere at equatorial plane. This corresponds to Region 4 in Fig.Q 

2.4.1 Density and velocity 



The disc density distribution is given by 

1 -(^4-^ 



p d (w, z) = E (w) 



2-Kh(w) 



(7) 



where w — \/ x' 2 + y 2 , h, z and E are the distance from the sym- 
metry axis, the scale height, the distance from the disc plane, and 
the surface density at the mid-plane, respectively. The mid-plane 
surface density and the scale height are given as: 

5M d 



EH 



-3/4 



(8) 



where i?do and Aid are the disc radius and the disc mass respec- 
tively, and 



h(w) = 0.05i?doW 9/8 . 



(9) 



With these parameters, the disc is slightly flared. The inner radius 
of the disc is set to Rdi = Rmi ■ The disc mass, A'ld, is assumed to 
be 1/100 of the stellar mass (M»), and the outer disc radius (_Rdo) 
is 100 au. The velocity of the gas/dust in the disc is assumed to be 
Keplerian. 

Because of the geometrical constraints imposed by the magne- 
tosphere and the disc wind, the finite height of the a disc causes an 
undesirable interface problem at the boundary; hence, the inner part 
(Rdi < w < -Rwo) of the disc is replaced by the geometrically thin 
(but opaque) disc. Since the emissivity of the disc at optical wave- 
lengths is negligible, and the disc scale height is relatively small 
in the inner part of disc, this simplification is reasonable: the most 
important factor in determining the shape of the Ha profiles is the 
finite height of the disc at large radii for large inclinations. 

2.4.2 Dust model 

In order to calculate the dust scattering and absorptio n cross s ection 
as a function of wavelength, the optical constants of iDraine & Led 
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Figure 2. The dependency of the disc-wind density and velocity on the wind acceleration parameter j3. The wind density p (right panel) and the poloidal 
velocity component v p (left panel) along the stream line starting from the mid point of the wind launching zone, i.e. (w, z) = («> m id> 0) where M) m jd = 
(^wi + Rwo) /2, are shown as a function of the distance (I) from the wind launching point (c.f. equations l5landltjl. The azimuthal velocity component (w^), 
which is independent of (3 (c.f. equationl4l. is also shown in the left panel for comparison. The density is normalised to the density po at the wind launching 
point for the /3 = 1.0 case. The v p reaches the terminal velocity by 100 au for all /3. In the far field (I > 10 au), the density is approximately proportional to 

~r 2 . 



1 1984) for amorphous carbon grains and Hannei 1 1988) for silicate 
grains are used. T he model uses the "large grain" dust model of 
I Wood etHil 120021) in which the dust grain size distribution is de- 
scribed by the following function: 



n (a) da — (Cc + Csi) a p exp 



da 



(10) 



where a is the grain size restricted between a m i n and a max , and 
Cc and Csi are the terms set by requiring the grains to com- 
pletely deplete a solar abundance of carbon and silicon. The pa- 
rameters adopted in our model are: C c = 1.32 x 1(T 17 , C Si = 
1.05 x 10" 17 , p = 3.0, q = 0.6, auxin = 0.1 ^m, a max = 
1000 fim, and a c = 50 urn. This corre sponds to Model 1 of 
the dust model used by IWood et alj 120021) . See also their Fig. 3. 
The relative numbers of each grai n is assumed to be that o f so- 
lar abundance, C/H~ 3.5 x 10~ 4 lAnders & Grevessdll989h and 
Si/H~ 3.6 x 10" 5 <Grevesse&Noelslll993l) which are similar 
to values found in the ISM model of Mathis, Rumol & Nordsieck 
Jl977l) and lKim. Martin & Hendrvl Jl994h . 



3 THE RADIATIVE TRANSFER MODEL 



We have extended the TORUS radiative transfer code iHarriesl 
l2000HKurosawa et all2004tlSvmington et all2005al) to include the 
multiple circumstellar components described above. In previous 
calculations (Svmington et al. 2005a), the model was used with a 
three-dimensional (3-D) adaptive mesh refinement (AMR) grid to 
investigate the line variability associated with rotational modula- 
tion of complex geometrical configurations of magnetospheric in- 
flow (see also iKurosawa. Harries & Svmingtonll2005l) . We modi- 
fied the code to handle a two-dimensional (2-D) density distribu- 
tion, and restricted our models to be axi-symmetric. Note that the 



velocity field is still in 3-D - the third component can be calculated 
by using symmetry for a given value of azimuthal angle. 

The computation of Ha is divided in two parts: (1) the 
source function calculation (S v ) and (2) the observed flux/profile 
calculation. In th e first process, w e have used the m ethod 
of iKlein & Castod <1978l) (see also iRvbicki & Hummed Il978t 
lH^rtmMiriT^1^994l) in which the Sobolev approximation method 
is applied. The population of the bound states of hydrogen are as- 
sumed to be in statistical equilibrium, and the gas to be in radiative 
equilibrium. Our hydrogen atom mod el consist of 1 4 bound states 
and a continuum. Readers are refer to lHarrieslfcOOOl) for details. 

Monte Carlo radiative transfer (e.g. under the 

Sobolev approximation, can be used when (1) a large velocity gra- 
dient is present in the gas flow, and (2) the intrinsic line width is 
negligible compared to the Doppler broadening of the line. In our 
earlier models lHarriej2000l : ISvmington et all20055) . this method 
was adopted since these conditions are satisfied. However, as noted 
and demonstrated by Muzerolle et al. 1 2001), even with a moderate 
mass-accretion rate ( 10~ 7 Mq yr -1 ), Stark broaden ing becomes 
important in the optically thick Ha line. iMuzerolle et all (2001) 
also pointed out that the observed Ha profiles from CTTS typi- 
cally have the wings exten ding to 500 km s -1 (e.g. [Edwards et alj 
ll994tlReipurth et alll99f3) which cannot be explained by the infall 
velocity alone. 

We have implemented the broadening mechanism following 
the formalism described bv lMuzerolle et all fcOOll) . First, the emis- 
sion and absorption profiles are replaced from the Doppler to the 
Voigt profile, which is defined as: 



H(a,y) 



(y - y') + a ' 



■dy 



(11) 



where a = r/47r Ai/p, y = ( v — up) /Az/p, and y — 
(V — vo) /A^d (c.f. Mihalas 1978). uo is the line centre fre- 
quency, and Ai^d is the Doppler line width of hydrogen atom (due 
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to its thermal motion) which is given by Ai/d = (2kT /mu) 1 ^ 2 X 
(fo/c) where run is the mass of a hydrogen atom. The damping 
constant T, which depends on the physical condition of the gas, is 
parameterised by Vernazza, Avrett & Loeser 1 1973) as follows: 



nm 



10 16 < 



T 



5000 K 



2/3 



(12) 



where nui and n e are the number density of neutral hydrogen 
atoms and that of free electrons. Also, C ra d, CVdw an d Cstark 
are natural broadening, van der Waals broadening, and linear Stark 
broadening constants respectively. We adopt this parameterisa- 
tion along with the values of broadening constants for Ha from 
iLuttermoser & Johnsoijl 19921) . i.e. C ra d = 6.5x 10~ 4 A, C v dw = 
4.4 x 10" 4 A and C Sta rk = 1-17 x 10" 3 A. In terms of level 
populations and the Voigt profile, the line opacity for the transition 
i — > j can be written as: 



XI 



- fij rij 



H(a,y) 



(13) 



where fij, rn, rij, (j^and g 3 are the oscillator strength, the popula- 
tion of i-th level, the population of j-th level, the degeneracy of the 
i-th level, and the degeneracy of the j-th level respectively. m e and 
e are the electron mass and charge (c.f. Mihalas 1978). 

We further modified TORUS by replacing the Monte Carlo line 
transf er algorithm with a direct integration method (c.f. [Mihalas 
Il978h for computing the observed flux as a function of frequency. 
The integration of the flux is performed in the cylindrical coor- 
dinate system (p, q, i) which is obtained by rotating the original 
stellar coordinate system (p, (j), z) around the y axis by the incli- 
nation of the line of sight. Note that the t-axis coincides with the 
line of sight with this rotation. The observed flux (F„) is given by: 



F v = 



4ird 2 



p sin q I v dq dp 



(14) 



where d, p max , and J„ are the distance to an observer, the maxi- 
mum extent to the model space in the projected (rotated) plane, and 
the specific intensity (/„) in the direction on observer at the outer 
boundary. For a given ray along t, the specific intensity is given by: 



I v = Joe To ° + 



/ °° S„ (t) e~ T dr 
Jt 



(15) 



where 7q and S„ are the intensity at the boundary on the opposite 
to the observer and the source function ("q^/xv) of the stellar atmo- 
sphere/wind at a frequency v. For a ray which intersects with the 
stellar c ore, Jp is computed from the stellar atmosphere model of 
iKuruczl 11979b as described in Section l2~Tl and lo = otherwise. 
If the the ray intersects with the hot ring on the stellar surface cre- 
ated by the accretion stream, we set 7o = B v (T r i ng ) where B v 
is the Planck function and T r i ng is the temperature of the hot ring. 
The initial position of each ray is assigned to be at the centre of the 
surface element (dA — p singdqdp). The code execution time is 
proportional to where n p and n q are the number of cylin- 

drical radial and angular points for the flux integration, and n v is 
the number of frequency points. In the models presented in the fol- 
lowing section, n v = 180, n q — 100, and n„ = 101 are used 
unless specified otherwise. A linearly spaced radial grid is used for 
the area where the ray intersects with magnetosphere, and a loga- 
rithmically spaced grid is used for the wind and the accretion disc 
regions. 

The optical depth r in eQuation fT5l is defined as: 



r(t) 



where x„ is the opacity of media the ray passes through, is 
the total optical depth measured from the initial ray point to the 
observer (or to the outer boundary closer to the observer). Initially, 
the optical depth segments dr are computed at the intersections of a 
ray with the original AMR grid in which the opacity and emissivity 
information are stored. For high optical depth models, additional 
points are inserted between the original points along the ray, and r\ v 
are Xv values are interpolated to those points to ensure dr < 0.05 
for the all ray segments. 

For a point in the magnetosphere and the wind flows, the emis- 
sivity and the opacity of the media are given as: 



Xv = Xc + xf + °c 



(16) 



where 77^ and n] 1 are the continuum and line emissivity of hydro- 
gen, x!?, Xl > an( l °cs are the continuum, line opacity (equation ll3t 
of hydrogen, and the electron scattering opacity. Similarly, for a 
point in the accretion disc, 



'<)<■ 



= 



dust 1 dust 
"'abs ~r ™sca 



(17) 



where KaJJs* and k^ 1 " a st are the dust absorption, and scattering opac- 
ity which are calculated using the dust property described in Sec- 
tion l2.4.2l We neglect the emissivity from the disc: since the disc 
masses of CTTS are rather small (~ 0.01 Mq) and low temper- 
ature (< 1600 K), the continuum flux contribu tion at Ha wave - 
length is expected to be neglig ible (e.g. lChiang & Goldreicrl 19971) . 



4 RESULTS 

4.1 Magnetosphere only models 

Using the reference model parameters (Table0 for the central star 
and the magnetosphere, we examine the dependency of Ha on the 
temperature (T max ) of the accretion flow and the mass accretion 
rate (M acc ), as done bvlMuzerolle et ail 1200 lh for H/3. The tem- 
perature structure of Hartmann et al. 1 1994) is adopted here, and 
the parameter T max sets the maximum temperature of the accre- 
tion streams. We note that here we have c omputed the hot ring 
temperature self-consistently, whereas Muzerolle et al. (2001) used 
a constant hot ring temperature (8000 K) for most of their mod- 
els. The accretion luminosity (L acc ) for models with A7 aC c = 
10~ r Mq yr _1 is about a half of the photospheric luminosity. 

Ha line profiles for a range of T max and A7 acc are presented 
in Fig.0 The overall dependency on T max and M acc is similar to 
that of Muzeroll e et alj 1200 lh . In general, the line strength weak- 
ens as the accretion rate and the temperature drop. The red-shifted 
absorption becomes less visible for the higher accretion rate and 
higher temperature models as it becomes filled-in by the stronger 
Stark-broadened line-wing emission. 

We demonstrate this effect in Fig. which shows an ex- 
ample for an Ha model with T max = 7500 K and M aC c = 
10 _7 MQyr -1 , both with and without damping. Although the 
maximum flux of the model with broadening is almost identical 
to that of the model with no damping constant (F = 0), a signif- 
icant increase of the line flux in both red and blue wings is seen. 
The weak red-shifted absorption component (which is a signature 
of the magnetospheric accretion) is weakened or eliminated by the 
flux in broadened wing. 
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Figure 3. Ha model profiles for wide ranges of mass accretion rate (M acc ) 
and temperature (T max ). The profiles are computed using only the magne- 
tospheric accretion flow (i.e. no outflow). All the profiles are computed us- 
ing the parameters of the reference model (TablelT] and inclination i = 55° . 
The temperature (indicated along the vertical axis) of the model increases 
from top to bottom, and the mass accretion rate (indicated by the values in 
Mq yr~ 1 along the top) decreases from left to right. 
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Figure 4. The effect of line broadening on Ha. A model computed with a 
damping constant (r), described in Section [3] (solid), is compared with a 
model with no damping, T = (dashed). Both models are computed with 
Tmax — 7500 K, i = 55°, and the reference parameters given in TableQ 
The two models have similar peak flux levels (around V ~ km s — 1 ), 
but the total flux and the EW of the line increased drastically for the model 
with the damping constant. The broad wings extend to ~ ±800 km s — 1 . 
The redshifted absorption feature seen (very weakly) in the T = model is 
not seen in the model with the broadening. 
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Figure 5. The dependency of the Ha profile on inclination (i). The profiles 
are computed with a magnetospheric accretion flow only using the reference 
parameters given in Table HI and T max = 7500 K. The solid, dotted, and 
dashed lines are for i = 10°, 55°, and 80° respectively. As the inclination 
becomes larger, the peak flux and the equivalent width of the line become 
smaller. 

Table [2] shows the equivalent width (EW) for the models in 
Fig. The EWs for half the models fa ll within the range o f 
EWs (~ 3 to ~ 160 A) measured by lAlencar & B asri (2000). 
The EWs of models with the lowe st mass accretion rates fall be- 
low the minimum EW observed by lAlencar & Basil j200C), while 
several models would be designated as weak-line d T Tauri stars 
(WTTS) using the traditio nal 10 A cut-off (e.g. iHerbig & Bel] 
1988). IWhite & Basri 120031) empirically showed that the full width 
of Ha at 10 per cent of its peak flux (10 per cent width) is a better 
indicator for an accretion than the EW criteria. They proposed that 
if a T Tauri star shows a 10 per cent width greater than 270 kms" 1 , 
the star should be a CTTS. The 10 per cent widths of the model 
H a profiles (Fig.|3j are summarised in Table|3] Using the criteria 
of IWhite & B asri (2003), half of the models shown in the figure 
can be classified as CTTS, and the other half as WTTS. How- 
ever, the profiles corresponding to the WTTS designations have 
an inverse P-Cygni m orphology that is rarely seen in observations 
iReipurth et alll99d and Sectionl5~TV 

The dependency of the line profile on inclination angle (z) 
is demonstrated in Fig. |5] The model has r max = 7500 K and 
Ma.cc = 10 -7 Mq yr _1 . The figure shows that the peak (nor- 
malised) flux decreases as the inclination angle increases, as does 
the EW. The line flux and the EWs decrease mainly because (1) the 
fraction of the accretion stream blocked by the photosphere in- 
creases as the inclination increases (c.f. Fig. Q lower-left part of 
the accretion stream i.e. z < and x < 0), and (2) the fraction of 
the accretion stream blocked by the inner part of the accretion disc 
increases as the inclination increases (c.f. Fig.Q lower-right part 
of accretion stream i.e. z < and x > 0). While the photosphere 
mostly blocks the blue-shifted part of the stream, the inner disc 
blocks the red-shifted part. The line EW in Fig.|5]changes from 32 
to 21 A as the inclination changes from 10° to 80°, i.e. a fractional 
change of 0.65. On the other hand, the 10 per cent width of the line 
changes from 290 to 470 km s _1 ; hence, the fractional change is 
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Table 2. Summary of Ha equivalent widths (A) from the magnetospheric 
accretion flow models shown in Fig. [3] Note that the equivalent widths are 
positive when the lines are in emission. 
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Table 3. Summary of 10 per cent widths (km s 1 ) from the magnetospheric 
accretion flow models shown in Fig.l3l 
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1.6. A similar EW dependence on inclination angle is found for the 
models with different magnetospheric temperatures and the mass- 
accretion rates. 

Because of the geometry of the magnetospheric accretion (c.f. 
Fig.0 and of the presence of the gas with the highest velocity close 
to the stellar surface, the highest blue-shifted line-of-sight velocity 
is seen only at the high inclination angles. For example, the ac- 
cretion stream close to the photosphere (~ 1_R» above the photo- 
sphere) on the upper-left quadrant in Fig.Q(s < and z > 0) 
will be more aligned with the line-of-sight of an observer in the 
upper-right quadrant in the same figure (x > and z > 0) as the 
inclination increases. This results in a slight increase in the blue 
wing flux for a higher inclination model as seen in Fig. [5] 

Our models show a blue-shifted line peak and a blueward 
asymmetry caused by occultation of the accretion flow by the stellar 
photosphere (see also Hartmann et al. 1994; Muzerolle et al. 2001). 
However, Alencar & Basri 2000 (see their fig. 9) found a substan- 
tial fraction of the observed Ha profil es show a red-shif ted peak. 
Furthermore, a recent study by Aooenzeller et al. (2005) demon- 
strated that the equivalent width of Ha from CTTS increases as 
the inclination angle increases. Clearly a magnetospheric accretion 
only model cannot explain the full properties of the Ha profiles in 
CTTS. 

4.2 Disc-wind only models 

In this section, we examine profiles produced using the disc wind 
outflow model described in Section 1231 The parameters used for 
the central star are as in Section l2~T1 and the disc-wind parameters 
are summarised in Table |4] Although the line is potentially sen- 
sitive to the temperature structure of the wind, determination of a 
self-consistent wind temperature is beyo nd the scope of this paper. 
Readers are refereed to Hartman n et alJ 11982T) in which the wind 
temperature structure is determined by balancing the radiative cool- 
ing rate (assuming optically thin) with the MHD wave heating rate. 
We pragmatically assume that the wind is isothermal at Twind- 
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Figure 7. The dependency of the Ha profile on the source displacement 
parameter d. The line profiles are computed with the wind parameter 
T wimJ = 9000K and /3 = 0.5 for i = 55°. The profiles in the figure 
use d = 2 R, (dashed), 7 R, (dash-dot) and 22 i?* (solid). All the other 
parameters are same as in Fig. |^] As the wind becomes more collimated 
(larger d), the wind emission becomes stronger. The P-Cygni absorption 
disappears for the model with a small collimation (e.g. d = 2 R* case) 
since the opening angle of the disc wind is larger than the inclination angle. 



Initially we examine the characteristics of the Ha profile as 
a function of the wind acceleration parameter (3 and the isother- 
mal wind temperature T w i n d. The mass-loss rate is kept constant at 
Mwmd = 10 -8 Mq yr -1 . Fig.|5|shows the model profiles com- 
puted for the wind temperatures between 7000 K and 9000 K, and 
those with j3 between 0.2 and 2.0, computed at inclination i = 55° . 
The morphology of the profile exhibited by the model changes from 
FU-Ori type (c.f. lReipurth et al!T9 96^ for the models with small (3 
to the stronger emission ty pes, which are m ore commonly seen in 
the observed Ha (c.f. Reiourth et al. 1996), as the value of (3 in- 
creases. Classical P Cygni profiles are prominent in the models with 
lower values of T w i n( j and (3 (a colder and faster accelerating wind). 
As the value of (3 increases, the position of the blue-shifted absorp- 
tion component moves toward the line centre. The blue edge of the 
absorption component in the fastest accelerating wind (J3 = 0.2) 
model is almost at the terminal velocity {vac ~ 160 km s" 1 ) of the 
disc wind (c.f. Fig.|2j, that for the slowest wind acceleration model 
(J3 — 2.0) is located close the line centre (~ 20 km s _1 ). 

The intensity of the emission component increases as j3 in- 
creases for a fixed T w j n a. This is mainly because the density near 
the photosphere increases as the value of (3 increases (Fig. |2j; 
hence, the emissivity increases. Most of the profiles show a red- 
ward asymmetry, mainly due to the presence of the classical P- 
Cygni absorption. However we find that the disc-wind models with 
the fastest acceleration (/3 = 0.2 and 0.5) and viewed at low in- 
clinations, show a strong blue asymmetry due to occultation of the 
receding half of the outflow by the optically thick disc. 

We now examine the sensitivity of the model profiles to the 
wind collimation or the source displacement parameter d. As we 
can see from Fig.Q the disc wind geometry/collimation is deter- 
mined by d and the wind launching radii (i? w i and R wo )- Here, 
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Figure 6. Ha profiles computed with the disc-wind only. The wind mass-loss rate is fixed at Af w ; n( j = 10 -8 Mq yr — 1 , and the line profiles are computed 
with different combinations of the wind acceleration rate (j3) and isothermal disc-wind temperature (T w ; n( j). The wind emission grows as the values of /3 and 
T w ; n( j increase. The position of the wind absorption moves toward the line centre as j3 becomes larger. 



we fix the values of 7? w i and _R wo as in Table [4] and inspect 
the model dependency on d alone. The minimum wind opening 
angle (measured from the symmetry axis) is given by 8 opcri = 
arctan (R„i/d). Fig.|7|shows the profiles computed for d = 2 _R*, 
7 R, and 22 R* for the inclination angle i = 55°. As the wind be- 
comes more collimated (larger d), the wind opening angle becomes 
smaller and the wind density is enhanced towards the symmetry 
axis (z); hence, the wind emission increases. When the collimation 
of the wind is small (as in the d — 2 i?* case), the P-Cygni ab- 
sorption is absent from the profile because the opening angle of the 
disc wind is greater than the inclination angle. For the disc wind 
models presented in Fig.[6| we use d — 22 R, in order to have 
a moderate amount of wind absorption at a moderate inclination 
angle (i ~ 60°). 



Blandford & Pavnd JT982T) showed that in order to launch a 
wind from a disc plane magneto-centrifugally, the launching angle 
(6l, measured from the equatorial disc plane) must be less than 
60°. In our models, the wind launching zone is restricted in be- 
tween R„i < Wi < Rwo (see Fig.Qand Table0 with _R w i = 3 /?* 
and _R wo = 1 au. This sets the upper limit of the source dis- 
placement distance (d) to be around 5.2 R, if we want to satisfy 
the launching angle condition of Blandford & Pavne 1 1982) for all 
the stream lines in our disc wind model. In the models presented 



in Fig. [6] (and later in Fig. [5] in next section), we have adopted 
d = 22 _R* instead. This means that some part of the disc wind 
(3-R* < Wi < 12.7 i?*) may not be able to be launc hed f r om th e 
disc plane according to the model of Blandford & Payne] ll982t) . 
We have used a larger launching angle or a larger source displace- 
ment distance for the following two reasons. 

First, we want to simulate a more collimated wind (toward 
the symmetry axis). Since the stream lines in our simple disc wind 
are represented by straight lines, the wind never bends toward the 
symmetry axis (z-axis) as the distance from the disc increases un- 
like the model of Blandford & Pavne 1 1982). We therefore used a 
larger launching angle to simulate the stream lines in the far field 
(at a larger distance from the disc plane) in their model. This is an 
intrinsic problem with approximating the stream lines with straight 
lines. 

Second, the larger source displacement model would provide 
a greater consistency with the fraction of Ha profiles affected by 
wind absor ption seen in ob servations. Accordin g to the o bserva- 
tions of Reipurth et alJ ll99o) and Alencar & Basri 1 2000), more 
than 50 per cent of CTTS show wind absorption features in their 
Ha profiles. In our model, if the launching angle is too small 
(#l < 60°), the wind absorption feature would not be seen or 
would be very small for profiles computed with i < 60°, as 
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Figure 8. The effect of the wind mass-loss rate. The profiles computed 
with (T wind , f3) = (9000K, 0.5), and M wind = 2.0 X 10" 8 M yr" 1 
(dashed), 1.0 xlO~ 8 M Q yr- 1 (solid), and 0.5 X 10~ 8 Mq yr -1 (dash- 
dot) are shown above. With all other parameter fixed, the mass-loss rate 
behaves as a scaling factor for the wind density (c.f. equationlol. 



Table 4. A summary of the reference disc-wind model parameters used in 
Sections|^|and|Sl See also SectionIO 
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demonstrated in Fig.Q The situation will however change if the 
field/stream line is not straight and bending towards the symmetry 
axis as the wind travels higher above the disc. 

In addition to the models given in Figs.[f>|and[9] we have com- 
puted the models with d = 5 R * (not shown here) which satisfy 
the wind launching condition of iBlandford & Pavnd |l982) in all 
parts of the disc wind zone in our model. These models showed 
that the wind absorption and emission features weaken, but over- 
all morphology dependence of the profiles on the wind temperature 
and the wind acceleration rate is similar to what we found in the 
models shown in Figs.[6|and[9| The minimum inclination angle at 
which the wind absorption becomes important is larger than that 
of the models with a larger source displacement (d = 22 R,) as 
expected. 

Finally, the effect of varying the wind mass-loss rate is shown 
in Fig. [8] With all other parameters fixed, M w i n d acts as a scaling 
factor for the wind density (equation |6j. The most important pro- 
cess that populates the n = 3 level is recombination, and therefore 
the line emissivity varies as the square of the density (assuming that 
the ionisation fraction remains constant); hence the line is very sen- 
sitive to changes in Af w i n( j. We indeed find that the line flux of the 
models is approximately proportional to the square of the mass-loss 
rate. 



4.3 Disc-wind-magnetosphere hybrid models 

We now consider models computed with a combination of the disc- 
wind and magnetospheric accretion components. The parameters 
used for the magnetosphere are same as in Fig.EKwith i = 55° ), 
and the mass-loss rate of the disc wind is M w i n d = 10~ 8 Mq yr -1 
(i.e. /i = M w i n d/M a cc = 0.1). Fig. [5] shows the models profiles 
computed (at i — 55°) using the same ranges of the disc-wind tem- 
perature and wind acceleration parameter as in the previous cases. 

As in the disc-wind only models (Fig.0, the location of the 
absorption component moves toward the line centre as the value 
of (3 increases. The figure shows that in the models with smaller 
(3 and T w i n d, the line emission from the magnetosphere domi- 
nates, while this situation reverses for the models with larger j3 
and T wind . The model with (T wind , (3) - (9000 K, 2.0) results 
in an Ha profile that i s far too strong to be c ompatible with obser- 
vations l Reipu rfhet all 19961 lAlencar & B asri 2000 ) ; however, the 
lower mass-loss rate models (while keeping fi constant) produce 
profiles with the line strengths comparable to observations (c.f. Ap- 
pendix^. Although the models in Fig.|5|are computed with limited 
ranges of T w i n d, /3 and a fixed i, the resulting profiles exhibit a wide 
variety of line profile shapes, many of which are simila r to the typ es 
of Ha profiles seen in the observations (c.f. Rei purth et alll 996). 

In order to quantify the relative contributions of the wind 
emission and the magnetosphere, we have computed the ratio of 
the EW of the hybrid model in Fig.|5|to the EW of the magneto- 
sphere only model (Table [5J- This ratio falls below 1.0 when the 
wind only contributes to the line absorption, but not to the emis- 
sion. In some models (e.g. with T w i n d = 8000 K and (3 = 2.0), 
the contribution of the wind to the EW is much larger than that of 
the magnetosphere (EW ratio is 3.2). This clearly demonstrates that 
the difficulty of using the EW of Ha alone as an accretion measure. 
We have also seen in Fig.|6|that disc wind models without a mag- 
netosphere can produce an Ha line with significantly large EW. 

An examination of the profiles in Figs.|6|and|9|reveals some 
problems with uniqueness that must inevitably occur when try- 
ing to assign an individual profile as wind or accretion dominated 
i.e. very similar profiles may arise from very different circum- 
stellar geometries. For example, the disc-wind only model with 
(Twind, P) = (8000 K, 1.0) in Fig.|6]and the accretion dominated 
model with (T wind , P) = (9000 K, 0.5) in Fig.^gfve very similar 
profiles. 

The effect of changing the ratio of the wind mass-loss rate 
to mass-accretion rate (fi = M w i nd /M acc ) is demonstrated in 
Fig. 1101 With a fixed value of mass-accretion rate (M acc = 
10 _7 Mq yr _1 ),the mass-loss rate is varied. The figure shows that 
as n increases the P-Cygni absorption deepens only slightly, and 
the position of minimum flux in the absorption trough appears to 
remain at same location. On the other hand, the flux peak increases 
and the line wings bec ome stronger as a incr eases. Although MHD 
models suggest (e.g. Koniel & Pudritz 2000) that \i ~ 0.1, the fig- 
ure suggest that models with a relatively wide range of \i (e.g. be- 
tween 0.05 and 0.2 in this case) produce profiles similar to obser- 
vations. 

By inspecting the model Ha images computed (Fig. II It . we 
have found that the extent of the Ha emission region is relatively 
small (i.e. h ~ 0.1 au or 10 R* where h is the vertical distance 
from th e disc) when c ompared to the spectro-astrometric observa- 
tions of lTakami et all 1200 3l) for RU Lup and CS Cha which show 
1-5 au scale outflows 1 . In our model, the wind emission flux drops 



Their observations also show some objects (e.g Z CMa and AS 353A) 
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Figure 9. The disc-wind-magnetosphere hybrid model. These models (solid lines) have the same parameters as Fig.[5| but also include the magnetospheric 
accretion flow (T ma x = 7500 K and M aC c = 10 -7 Mq yr -1 ). While the emission from the magnetosphere dominates for the models with smaller /3 and 
T w ; nc [, the emission from the wind dominates the profiles for models with larger /3 and T w ; n( j. The contribution of the magnetosphere can be judged from the 
profiles computed with the magnetosphere only (dotted lines). 



by ~ 10 4 by the time the wind reaches h w 0.1 au. The extent 
of the line emission region becomes slightly larger (~ 0.4 au) if 
a slower wind acceleration rate (e,g, f) = 4) is used. The discrep- 
ancy between the observations and the model is perhaps caused by 
the lack of the material connecting the disc wind to a collimated jet 
in our model, which may be excited by alternative heating mecha- 
nism (e.g. shocks) thus producing observable Ha emission at large 
distances from the central object. 



Table 5. The ratio of the line equivalent width (EW) of the disc-wind mag- 
netosphere hybrid models in Fig.|9]to the EW of the magnetosphere only 
model. 
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5 DISCUSSION 

5.1 The morphological classification scheme proposed by 
Reipurth et. al (1996) 

iReipurfh et alj fT996) proposed a two-dimensional classification of 
Ha emission profiles of T Tauri stars and Herbig Ae/Be stars. Their 



displaying the outflows in larger scale (> 50 au); howe ver, this could 
be fo rmed in shocks rather than MHD-wave heating (e.g. Hartmann et al. 
Il982l) or X-ray heating (e.g. lShang et all2002t) . 



classification scheme contains four classes (I, II, III and IV) differ- 
entiated by the ratio of the secondary-to-primary emission compo- 
nents in the profiles. Each class is divided into two sub-classes (B 
and R) which depends whether the secondary peak is on the blue 
or red side of the primary peak. Readers are referred to Fig. 4 of 
their paper. Fig. ll2l shows sample model profiles which are classi- 
fied according to the definition of lReipurth et alj 1 1996). The com- 
bination of the disc wind, magnetospheric accretion flow, and the 
accretion disc can reproduce all the classes of the profiles seen 
in observations. The corresponding parameters of the disc-wind- 
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Velocity (km s ) 

Figure 10. The effect of the mass-loss rate to mass-accretion rate ratio 
(p = Af w ; n( j/Af acc ). The Ha profile computed with only the magneto- 
spheric accretion (thick solid) is compared to the disc-wind-magnetosphere 
hybrid models for fj, = 0.05 (dash-dot), 0. 1 (solid) and 0.2 (dash). The mag- 
netosphere used here has T max = 7500 K and Af acc = 10 -7 Mq yr _1 . 
The temperature and the acceleration parameter of the disc wind are 
T\vmd = 9000 K and /3 = 0.5 respectively. All the models are computed 
with the inclination angle i = 55°. As fi increases the emission component 
becomes stronger and the wings becomes broader, as expected. The absorp- 
tion component is relatively insensitive to the value of fj, in this example. 
The wide range of fi reasonably reproduces the profiles similar to ones seen 
in observations (e.g. iReipurth et all 19961) . 
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Figure 11. An Ha image of the disc- wind hybrid model viewed at an incli- 
nation of 80°. The false-colour scale shows logarithmic (base 10) surface 
brightness (arbitrary units) and the axes are linear dimensions in au. The 
parameters used for the magnetosphere are as the reference model in Ta- 
ble[3 The wind temperature of T m ax = 7000 K and the wind acceleration 
parameter /3 = 1.0 are adopted here. 



Figure 12. Sample Ha model profiles which characterise the classification 
scheme of Reirjurth et al. 1 1996). The combination of magnetospheric ac- 
cretion flow, the accretion disc, and the disc-wind can reproduce the wide 
range of Ha profiles seen in observations. The model parameters used for 
each profile along with short comments are given in Table|S| The horizontal 
axes are velocities in km s , and the vertical axes are normalised flux. 



magnetosphere hybrid model used to reproduce the profiles in the 
figure are summarised in Table [5] along with brief comments on 
possible physical conditions which lead to the profiles in each 
class. Below we discuss the Reipurth scheme in the context of 
our models parameters. Readers are referred to Appendix lAl for 
the complete model profiles used in this discussion. Although in 
general Ha in the spectra of CTTS shows a significant amount of 
variablity, which may in some cases be ind uced by the rot ational 
motion of an inclin ed magnetosphere (e.g. I Johns & Basrj[l995l 
iBouvier et alJl2003h . we restrict our discussion in this section to 
the axi-symmetri c case only. 

Type I: Reipurth et al. ( 1996) found that Type I profiles (sym- 
metric around line centres) constituted 26 per cent of the 43 CTTS 
Ha profiles in their sample, making it the second most common 
morphological type. We find that Type I profiles are dominated by 
magnetospheric emission at rates in excess of ~ 10~ 8 Mq yr _1 . 
Type I profiles appear at a wide range of inclinations (Fig. [5}- An 
object that displays a Type I profile may have relatively weak disc 
wind, but in this case it must be viewed from the pole (in order that 
no wind material is seen projected against the stellar photosphere). 
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Type II-B: These profiles have a secondary blue peak that is 
in excess of half the strength of the primary peak, and comprised 
16 per cent of the Reipurth sample. We find that Type II-B profiles 
only occur at medium to high inclinations, and for systems with 
mass accretion rates of w 1CP 7 Mq yr _1 . If the system is viewed 
at high inclination, a Type II-B profile implies a fast wind acceler- 
ation, and conversely a moderate inclination angle suggests a more 
modest wind acceleration. 

Type II-R: These profiles are characterised by a secondary 
red peak that is in excess of half the strength of the primary peak. 
They are approximately as common as the Type II-B lines. For 
the M acc = 10~~ 7 Mq yr -1 a switch from Type II-B profiles to 
Type II-R occurs when the inclination angle changes from 55° to 
80°. This profile type is also seen at lower mass-accretion rates, but 
typically only in the models with the slowest accelerating winds. It 
is perhaps unsurprising that the Type II-B and Type II-R profiles 
are equally common, since (in general) the split between the profile 
types results from a geometrical (viewing angle) effect rather than 
a marked difference in physical parameters (mass accretion rate, 
wind acceleration etc). 

Type III-B: Profiles with a blue secondary peak that is less 
than half the strength of the primary peak are classified as Type Hi- 
fi . This type of profile is the most common in the Reipurth sample, 
comprising 14 out of 43 (33 per cent) CTTS. We find Type III-B 
profiles only in our Af acc = 10~ 7 Mq yr" 1 models. A fast wind 
acceleration and a moderate inclination are necessary, or a slower 
accelerating wind viewed pole-on. We do not see this profile type 
in our high inclination models. Although Type III-B profiles do not 
occur at Macc = 1-0 x 10 _8 MQyr _1 , we have examined the 
intermediate case of M acc = 3 x 10" 8 Mq yr" 1 , which may rep- 
resent a typical T Tauri mass accretion rate iGullbring et alll998t 
ICalvet et aj2004 . We find that the models with fast accelerating, 
high temperature winds do display Type III-B profiles. It is clear 
that our model, as it stands, requires above-average (but perhaps 
not atypical) mass-accretion rates in order to produce Type III-B 
profil es. It is less clea r whether the 14 Type III-B objects found by 
iReipurfh et alj fl996) constitute the highest mass-accretion rate ob- 
jects in their sample: the average EW for T Tauri stars with Type III- 
B profiles is 64 A, which is in the top quartile of their complete 
T Tauri star sample, perhaps indicating that Type III-B profiles are 
associated with high EWs, and therefore strong mass-accretion. 

Type III-R: This profile type has a red secondary peak that 
is less than half the strength of the primary peak, and only one 
object with this type of profile was observed in the Reipurth sam- 
ple (SZ Cha). We are able to reproduce this profile type, but only 
using a very narrow range of parameter space. In fact there are 
no Type III-R profiles in our model grid, but we obtain one by 
'tweaking' the magnetospheric temperature, and using an inclina- 
tion (85°) sufficient that the disc starts to obscure the central Ha 
emission. For completeness, the final model parameters we used to 
achieve the profiles are: very high inclination, T max ~ 8300 K, and 
M a cc ~ 10 -8 Mgyr^ 1 . It is interesting that the most infrequently 
observed profile type was in fact the most difficult for us to repro- 
duce. It appears that this profile morphology requires some obscu- 
ration effects by the dust disc, and therefore a high inclination; this 
naturally explains the rarity of the profile in the observations. 

Type IV-B: This profile classification is the classical P-Cygni, 
in which the blue-shifted absorption component has a sufficient ve- 
locity to be present beyond the emission line wing. It occurs in- 
frequently in observations (7 per cent of the Reipurth sample). No 
Type IV-B profiles were produced in our model grid, since the blue 
wing of the magnetosphere has generally such a large extent that it 



always exceeds the projected wind terminal velocity, and thus the 
bluemost extent of the absorption component. Only by running a 
model with a larger terminal velocity or with a smaller magneto- 
spheric temperature, were we able to achieve a classical P-Cygni 
shape. We suspect that the disc-wind geometry is not the best kine- 
matic model for producing this morphology. A more straightfor- 
ward explanation lies with a bipolar flow, viewed along its axis, or 
even a simple spherical wind. 

Type IV-R: This is the inverse P-Cygni profile type, and 
only occurs in about 5 per cent of Ha profiles. Type IV-R pro- 
files occur very infrequently in our hybrid model grid, with only 
the lowest mass-accretion rate models viewed at high inclination 
showing the characteristic shape (although using the EW criterion 
these models would not be classified as CTTS). The inverse P- 
Cygni profile requires that the observer view the hot accretion foot- 
prints through the fast moving magnetospheric accretion flow, and 
this means that the inclination angle is restricted to small range. 
We can reproduce Type IV-R profiles more readily if we adopt a 
cooler magnetospheric temperature than indicate d by the ac cretion- 
rate/temperature relationship presented bv lMuzerolle et all 1200 lh . 
The restriction in the inclination angle mentioned above would be 
relaxed if the axis of the magnetosphere was tilted with respect to 
the rotational axis since the line-of-sight from an observer to the 
hot rings will pass through different parts of the magnetosphere as 
the rotational phase changes. 

5.2 The inclination dependence of Ha profiles 

Recently, Appen zeller et alj 12005F) have shown that the EW of the 
observed Ha from CTTS is inclination dependent, with the EW 
increasing as the inclination angle increases. Next, we will exam- 
ine whether the hybrid model is consistent with the inclination de- 
pendency found in the observation. The left-hand panel in Fig. 1131 
shows the profiles computed with the disc-wind, magnetosphere 
hybrid model (Section f4.2> at inclination angles i — 10°, 55° and 
80°. The magnetospheric accretion component uses the parameters 
used as the reference model (Table 0. The figure shows that the 
absorption feature becomes stronger as the inclination increases, 
mainly because of the geometrical configuration. The optical depth 
to the observer becomes larger as the inclination angle increases, 
since the density of the disc wind increases towards the equatorial 
plane - as a result, the EW becomes smaller as the inclination angle 
increases: 41, 27, and 19 A for i — 10°, 55°, and 80° respectively. 
This is a robust result, and is insensitive to the adopted wind pa- 
rameters. Recall that this is the same EW/inclination dependency 
as the magnetosphere only models (Fig. [5}- 

Clearly the inclination dependen cy of the l ine strength is in 
contradiction with the Appenzeller et al. (2005) result, and thus 
merits further investigation. In the previous sections, we have con- 
sidered Ha line formation in a disc wind, which mimics the den- 
sity distrib ution of the mag neto-centrifugal launched jet models 
(e.g. lKrasnopolskv et aljE003T) . The density in this model is more 
concentrated toward the equatorial plane in the region where the 
line emissivity is highest, although further out the material becomes 
collimated into a bipolar jet, this would not contribute significantly 
to the line emission in our models. 

As an alternative, we consider a wind model in which the out- 
flow density increases towards the pole, rather than the equator. 
In this model, the disc-wind is replaced by a stellar wind, prop- 
agating only in the radial direction, accelerating under the classi- 
cal beta-velocity law (c.f. ICastor & Lamerdll979l) i.e. v r ( r) = 
Voo (1 — Rt/r) 13 where Voo and R, are the wind terminal velocity 
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Table 6. The summary of model parameters for the profiles in Fig. ll2l and brief comments. The temperatures are in 10 3 K. M acc and M wind are in Mq yr~ 1 . 



Class 


i 




Tin ax 


Af w ind 




/3 


Comment 


I 


55° 


10" 


-7 


7.5 








Accretion dominated. Wide range of inclination. 


II-B 


55° 


io- 


-7 


7.5 


IO" 8 


8.0 


1.0 


Wide range of wind acceleration rates. Mid-high inclination. 


II-R 


80° 


10" 


-7 


7.5 


IO" 8 


7.0 


2.0 


Slow wind acceleration rate. High inclination. 


III-B 


55° 


io- 


-7 


7.5 


IO" 8 


9.0 


0.5 


Fast wind acceleration rate. Mid inclination. 


III-R 


85° 


io- 


-8 


8.2 


IO" 9 


6.0 


1.0 


Mid wind acceleration rate. Very high inclination. 


IV-B 


30° 


io- 


-7 


5.5 


IO" 8 


6.0 


0.2 


Fast wind acceleration rate. Mid-high inclination. 


IV-R 


55° 


io- 


-9 


9.5 








Accretion dominated. Low mass-accretion rate. Mid inclination. 




Velocity (kms"') Velocity (km s ') 



Figure 13. The dependence of the Ha model profiles on the inclination angle. The profiles computed at i = 10° (solid line), 55° (dotted line) and 80° (dashed 
line) using (1) the disc-wind-magnetosphere hybrid model (left; c.f. Section l4"2l and (2) the disc, stellar wind, magnetosphere hybrid model (right panel; c.f. 
Fig. 1141 . The parameters used for the magnetosphere are same for the both models, i.e. T max = 7500 K and Af acc = 10~ 7 Mq yr~ 1 . The isothermal 
wind temperature and the wind acceleration parameter used in (1) are T w ; nt j = 7000 K and /3 = 1.0, and those in (2) are T w ; n( j = 9000 K and /3 = 1.0 
respectively. The wind mass-loss rate used in both models is M w ; n( j = 10 — 8 Mq yr _1 . All the other parameters are same as in Sections l4.3l and l4~2l 



and the stellar radius respectively. f3 is the wind acceleration pa- 
rameter similar to the one seen in the disc-wind model (c.f. equa- 
tion|5J- The wind emerges from the stellar surface, but restricted to 
be within ~ 30° from the symmetry axis hence forming the cone- 
shaped regions both above and below the poles. The density is com- 
puted from assuming a constant mass-loss rate per unit area over the 
stellar surface, requiring mass-flux conservation, and using the beta 
velocity law above. The schematic of the model configuration is 
shown in Fig.^| Again, we assume the wind temperature (T w i n d) 
is isothermal for simplicity. This model configuration is motivated 
by the recent study by M att & Pudritzl 120051) who demonstrated 
that a stellar wind along open magnetic field lines originating from 
the star may cause significant spin-down torque. 

The right-hand panel in Fig. ll3l shows Ha profiles computed 
using the alternative model at the same inclination angles used for 
the disc-wind-magnetosphere hybrid model. The mass-loss rate of 
the stellar wind in this model is 10 -8 Mq yr -1 , which is the same 
as the mass-loss rate used in the disc-wind magnetosphere hybrid 
model in the same figure. In both models, the mass-accretion rate 
of lO _7 M0yr _1 is used. The wind acceleration parameter (f3) 



and the wind temperature (T w i n d) used here are 1.0 and 9000 K 
respectively. All the other parameters are the same as in the disc- 
wind magnetosphere model. With this set of [3 and T w i nd , the wind 
emission is significantly larger than the emission from the magne- 
tosphere. Because of the geometry of the wind, the P-Cygni ab- 
sorption feature weakens as the inclination increases: the optical 
depth of the wind is much higher in the polar direction. For the 
same reason, the EW of the line increases as the inclination in- 
creases, as clearly seen in the figure. Specifically the values are 54, 
69, and 98 A for i = 10°, 55°, and 80° respectively. This inclina- 
tion dependency of EW is seen in the models with wide ranges of 
(3 and T w i n d, provided that the wind emission of the model is dom- 
inating the emission from magnetosphere. A similar inclination de- 
pendency of EW is found for lower wind mass-loss rate models 
(e.g. Mwind = 10~ 9 Mq yr -1 ) while keeping ji = 0.1. 

The A ppenzeller et alJi2005l) study used just 12 objects, and a 
significantly greater sample should be investigated before firm con- 
clusions are drawn. However our models indicate a significant dif- 
ference in the inclination dependence of the line morphology which 
may be simply tested against observation. For example, one possi- 
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Figure 14. The alternative outflow model. This model replaces the disc- 
wind with a stellar wind. The wind is launched from the stellar surface, and 
is limited to cones with opening angles of ~ 30° . The system consists of 
four components: (1) the continuum source located at the origin (o) of the 
cartesian coordinates (x, y, z) - the j/-axis is into the paper, (2) the magne- 
tospheric accretion flow, (3) the stellar wind outflow, and (4) the accretion 
disc. The density distribution is symmetric around the z-axis. 

ble way to distinguish the two different wind models is to examine 
the profile shape of a CTTS which is known to be viewed close to 
pole-on (e.g. a CTTS with a very low v sin i). If the profile contains 
a prominent blue-shifted absorption component (e.g. Type IV-B), 
the system is more likely to be associated with the bipolar stellar 
wind type of outflow, and not with the disc-wind. In reality, how- 
ever, the two types of winds may be present in the same obje ct, in a 
simila r fas hion to that seen in the models of iDrew. Proga & Stond 
( 1998) and lMatt & Pudritzl 120051 In these models, the fast stellar 
wind is present in the polar directions, while at the same time a 
slower, denser disc-wind is present near the equatorial plane. 



6 CONCLUSIONS 

We have presented disc-wind-magnetosphere hybrid radiative 
transfer models for classical T Tauri stars, and detailed studies of 
the Ha formation from their complex circumstellar environment. 
We found that the hybrid model can reproduce the wide variety of 
profile seen in observations (Figures I3ll6ll9landf!2l. 

Using the model res ults, we examined the Ha spectroscopic 
classification proposed by Reipurth et al. ( 1996), and discussed the 
basic physical conditions that reproduce the profiles in each classi- 
fied type (Section f5.1> . Using the different combinations of the in- 
clination (i), the mass-loss to mass-accretion rate ratio (/i) and the 
wind acceleration rate (/?), our radiative transfer model w as able to 
produce all 7 classes of profiles defined in Reicurth et al. 1 1996). 

One potential free parameter of our models is the ratio of 
mass-loss rate to mass-accretion rate (/x). M ost MHD simulations 
predict a value of n = 0.1 (c.f. Konigl & Pudritzl EboO). and we 
have adopted this figure for most of our models. The success of 
our profile in matching observations is therefore encouraging, and 
provides support for this canonical value of fi. We have computed 
sample models with other values of /i, (Fig. llOt and find that mod- 
els with 0.05 < fi < 0.2 are still consistent with observations. 
Wind-only or magnetosphere-only models may be representative 
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of individual line profiles, but neither can explain the complete 
Reipurth classification scheme - even magnetosphere dominated 
profiles (Type I) can be readily explained by the hybrid model with 
fi = 0.1. 

Unfortunately, the dependency of the line equivalent width on 
the inclination angle predicted by the disc-wind-magnetosphere hy- 
brid model does n ot agre e with the trend seen the observations of 
Appenzelle r et alj l2005t) . We may appeal to the small sample size 
of their study (12 objects) to question the validity of the EW effect, 
but we have also considered an alternative model in which the disc- 
wind is replaced by a biconical stellar wind, and have found that 
this alternative wind model can reproduce the line equivalent width 
dependency on the inclination CFig. 1131 . However, MHD simula- 
tions predict that bipolar or jet-like flows in CTTS are always ac- 
companied by a disc-wind originating near the star: the actual situ- 
ation therefore is likely to be a combination of the two flows. 

Although we have concentrated on CTTS here, the hybrid 
model is equally applicable to pre-main-sequence stars across 
the mass spectrum. Herbig Ae/Be star profiles are equally di- 
verse, and w hilst many of the profiles are indicative of a wind 
iFinkenzeller & Mundtl Tl984 some lines may be ex plained by 
magnetospheric accretion alone (e.g. Muzerolle et al. 2005). One 
possible advantage in applying the hybrid model to the Her- 
big Ae/Be systems is the growing observational data probing 
the line formation region and inner disc usin g techniques such 
as spect ro-astrometry (e.g . iBaines et alJl2005l) . sp ectropolarime- 
try (e .g. lvink et alll2005al) and interferometry (e.g. iMonnier et all 
2005). These high-resolution techniques may aid model fitting of 
individual objects by constraining some of the model parameters, 
in particular the inclination. 

Velocity broadened, asymmetric Ha profiles in brown dwarfs 
(BDs) are used as an indicator of ongoing magnetospheric ac- 
cretion, and the same models that are used for CTTS are ap- 
plied to the se systems in order to determine mass-acc retion rates 
(e.g. lLawson. Lvo & MuzerolldEooH : iMuzerolle et alJl2005l) . The 
line profiles however display the same wide variety of morpholo- 
gies as the CTTS and Herbig Ae/Be stars (although we note as yet 
no BD profile has been observed with red-shifted absorption). It is 
natural to assume that the same physics applies in these low mass 
objects, and that a hybrid model may be more successful in fitting 
these profiles, and will give greater insight into the circumstellar 
kinematics. 

There are several avenues of future work to consider. Here we 
have adopted a simple analytical form for the accretion and outflow 
velocity and density structure. A more self-consistent approach 
would be to use the results of MHD simulations as an input to the 
radiative-transfer modelling. The thermodynamics of the circum- 
stellar gas needs to be examined (c.f. Hart mann et alll994l :lMartin 
1996) and also computed self-consistently with the radiation field 
used in the radiative-transfer, a process that would decrease the free 
parameters of the model. We intend to calculate other observables 
from the model that may be compared with observation and should 
provide further constraints. The line polarisation signature, which 
arises from scattering of line and continuum photons in the circum- 
stellar dust, encodes both dynamical and geometric information, 
such as the size of the inner disk radius (e.g. IVink et alj l2005b: 
IVinketall2 005a). Spectro-astrometry may provide information on 
the system's inclinat ion, and the geometry o f the outflow on large 
scales (e.g. Takami et al. 2003: Baines et al. 2005). Interferometry 
also provides information on inner disc radii and inclinations, al- 
though currently only for the Herbig Ae/Be stars. 

We have concentrated on the formation of the Ha line, be- 
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cause it is the best studied CTTS emission line. However it is clear 
that simultaneous fitting of several emission line profiles should 
provide much stronger constraints. In particular it appears that the 
near-IR hydrogen lines (Pa/3 and Br7) are less affected by outflows 
than Ha. One could therefore use these lines to constrain the mass- 
accretion rate, and use Ha to obtain the mass-loss rate. A potential 
pitfall is the significant temporal variability that all CTTS show in 
their emission lines; the observational data on a particular object 
must be co-temporaneous. 
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APPENDIX A: Ha ATLAS OF THE 
DISC-WIND-MAGNETOSPHERE HYBRID MODEL 

Here we present our Ha model line profile grid computed us- 
ing the disc-wind-magnetosphere hybrid model in Section l4~3l All 
the models use the reference parameters for the magnetospheric 
geometry and the photosphere as given in Table \\\ in the main 
text. The ratio (n) of the mass-loss rate to mass-accretion rate 
is fixed to 0.1 for all the models. Three sets of model profiles, 
each with a different mass-accretion rate, are presented here: (1) 
Mace = lO^MQvr" 1 (Fig. ED, (2) M acc = lO^MQvr" 1 
(Fig. E3, and (3) A/ acc = lO^Mgyr" 1 (Fig. 1X51, The tem- 
peratures of the magnetospheres used for (1), (2) and (3) are 
7500 K, 8500 K and 9500 K respectively. The profiles are pre- 
sented for the inclination angles of 10°, 55°, and 80°. The 
disc-wind geometry used to compute these profiles are the same 
as in Section 14.31 These profiles are available electrically at 
|http : / /www .astro .ex.ac. uk/ people /rk/pr of iles| 



P=0.2 p=0.5 M-0 P=2-0 




V(kms"') V(kms"') V (tans' 1 ) V(tans"') 



P=0.2 p=0.5 p=1.0 P=2.0 




V(kms"') V(kms"') V(kms"') V (km s" 1 ) 



P=0.2 p=0.5 p=1.0 P=2.0 




-500 500 -500 500 -500 500 -500 500 

V(kms"') V(kms"') V (tans' 1 ) V(kms"') 



Figure Al. An Ha atlas of the disc-wind-magnetosphere hybrid model with 
M acc = lO _7 M0yi— 1 and M wind = 10~ 8 Mq yr _1 . Profiles are 
presented for inclinations of 10° (top figure), 55° (middle figure), and 80° 
(bottom figure). For each inclination, 12 profiles are given, for three differ- 
ent values of the wind temperature, and four different values of the wind 
acceleration parameter f3. The contribution of the magnetosphere can be 
judged from the profiles computed with the magnetosphere only (dotted 
lines). 



